Can a Large Neutron Excess Help Solve the Baryon Loading Problem 

in Gamma-Ray Burst Fireballs? 
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We point out that the baryon-loading problem in Gamma-Ray Burst (GRB) models can be ame- 
lioriated if a significant fraction of the baryons which inertially confine the fireball are converted to 
neutrons. A high neutron fraction in some circumstances can result in a reduced transfer of energy 
from relativistic light particles in the fireball to baryons. The energy needed to produce the required 
relativistic flow in the GRB is consequently reduced, in some cases by orders of magnitude. This 
could be relevant to GRB models because a high neutron-to-proton ratio has been calculated in 
neutron star-merger fireball environments. Significant neutron excess also could occur near compact 
objects with high neutrino fluxes. 



PACS numbers; 98.70.Rz, 14.60.Pq 

In this Letter we show how the baryon loading prob- 
lem can be alleviated in certain gamma-ray burst (GRB) 
models when significant numbers of baryons are con- 
verted to neutrons. Interestingly, many of the proposed 
GRB "central engines" involve compact objects which are 
themselves highly neutronized, or which are accompanied 
by intense neutrino fluxes. Weak interactions induced by 
these neutrino fluxes can result in significant proton-to- 
neutron conversion, especially if resonant neutrino flavor 
transformation takes place [|l]j2[|^. 

Inferences of the energetics and spectral observations of 
GRBs imply (i) total energies in gamma-rays approach- 
ing 10^'^ ergs for the most energetic events (in the absence 
of beaming) , and (ii) large Lorentz factors of the progen- 
itor fireball (7 ^ 10'^) (for a recent review, see Ref. [Q). 
Excessive baryon pollution of the fireball precludes at- 
tainment of these features for many GRB models. This 
is a consequence of the conversion of radiation energy in 
the electron/positron/photon fireball to kinetic energy 
in baryons [|[^. However, the relatively small cross sec- 
tions characterizing the interactions of neutrons with the 
electron/positron/photon plasma may afford a solution 
to this problem. 

This can be seen by considering the fictitious limit of 
completely noninteracting neutrons. Imagine that pro- 
tons inertially tether an electron/positron/photon fire- 
ball via photon Thomson drag on e^, which in turn in- 
fluences protons through Coulomb interactions. If these 
protons were suddenly converted to non-interacting "neu- 
trons" , then the fireball would expand relativistically, 
leaving behind the baryonic component. Real neutrons 
can approximate this limit as they interact with the elec- 
tron/positron/photon plasma only via the neutron mag- 
netic dipole moment. These cross sections are small com- 
pared to the Thomson cross section a^: neutron-electron 
(positron) scattering has a^e ~ 10~^gT H]; and neutron- 
photon scattering has aa^y ~ IO^^^ctt [HT 

However, the real limit on the efficacy of this mecha- 
nism is the strong interaction neutron-proton scattering 
which will dominate the energy transfer process when 
conversion of neutrons to protons is incomplete. There- 



fore, the degree to which the baryon loading burden can 
be lifted in our proposed mechanism will depend on the 
neutron excess in the fireball environment. Here we will 
measure the neutron content of the plasma in terms of 
the electron fraction Yg, the net number of electrons 
(rig- — ne+ ) per baryon, or in terms of the neutron-to- 
proton ratio Ye = l/(n/p+ 1). 

We note that although previous studies have invoked 
neutrino oscillations to attempt a baryon loading prob- 
lem solution P,p^, none has exploited the Ye-changing 
aspect of the weak interaction. 

To go beyond the simplistic picture of non-interacting 
neutrons, we can consider a two-component ((i) neutrons, 
and (ii) protons/e^ /photons) plasma in the context of 
a homogeneous fireball with initial radius, temperature, 
Lorentz factor, and electron fraction, i?o,21),7o & ^eo, 
respectively. Numerical and analytic work have shown 
the following simple scaling laws for such a configuration 
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For R < TjRo/-/a 
For R > i-jRo/lo 



7 = 70 R/Rf) 
T = To Ro/R 



7 = 77. 
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As a matter of convenience we will take the scaling to 
be such that 70 = 1. In Eq. (|^), the ratio of energy in 
radiation E to total baryon rest mass M is t] = E/M . 

One can relate rj to the entropy per baryon, s, using 
the number density of baryons N = Ph/rnp = Prad/mpT] 
where pt the baryon component rest mass energy density 
and nip is the proton rest mass. Using this relation, and 
noting that in terms of the proper entropy density S, the 
entropy-per-baryon is s = S/N, the relation between s, 
T], and the temperature Tq is s « 1250 r]{l MeV/To). 

For a large enough 77, baryon loading is unimportant 
g. In fact when ?7 > 10^{E/10^^ erg)i/3(10^ cm/i?o)^/^ 
the fireball becomes optically thin before transferring its 
energy to kinetic energy in baryons (here E is the total 
energy of the fireball). 

Written in terms of the time t as measured in a frame 
comoving with the fireball the above relations imply 
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T = To (2) 



for ^ < rj. Here the dynamic timescale is defined to be 
the initial hght crossing time, Tdyn = Ro/c. In fireballs 
resulting from neutron star mergers, for exam ple , Tdyn ^ 
2.7 X lO^^s, corresponding to an Rq of 8 km fl^ . 

A particle co-moving with the expanding plasma ex- 
periences a 4- acceleration a^^ with magnitude y/oFci^ — 
d'f/dR. As noted above, the force that drags the neu- 
trons along with the expanding plasma arises princiapally 
from n-p collisions. The relative contribution to the total 
force on the neutrons from collisions with electrons and 
positrons is roughly Fn-e/Fn-p ~ meneane/mpUpanp < 
10^^" (s/Ye) and is small for the conditions we con- 
sider. The neutron- photon cross section is small enough 
{an-y ~ 10-^^[E^/il MeV)]2 cm^ where is the pho- 
ton energy in the neutron rest frame [^) that n-7 inter- 
actions are negligible. 

The relations in Eqs. (|l|) imply that an inertial ob- 
server with time coordinate t' initially (at t' = 0) co- 
moving with the plasma sees the plasma accelerate ac- 
cording to 7f = t' /r^y^. Hereafter we adopt natural 
units where c = 1 . If we denote by t~^1^ the frequency of 
neutron/proton collisions (per neutron), we expect that 
the two components of the plasma will achieve a rela- 
tive velocity given by Wj-ci ~ 2rcoii/Tdyn, where the factor 
of 2 arises from the approximate angle independence of 
the neutron-proton scattering cross section and the near 
equality of the neutron and proton masses. An equiv- 
alent expression is found if one considers the force on 
the neutrons from collisions with protons . It is clear 
then that when Tdyn ^ TcoU the neutrons are coupled 
to the rest of the plasma. However, decoupling occurs 
as these two timescales become comparable. Since the 
baryon number density in the plasma frame decreases 
as e"^*/'^''*'", decoupling will occur quickly, i.e. on a 
timescale shorter than Tdyn. 

When significant decoupling occurs we can neglect the 
thermal contribution to the collision frequency and write 



'coll 



(9 X 10 



12 



7^3 

-l-j ^McV 



S5 



Vrcl YeCrio . 



(3) 



where S5 = s/lQ^ and crio is the neutron- proton cross- 
section in units of 10 fm^, and TmcV is the temperature 
in MeV. As the precise energy dependence of ctio is not 
important here, it suffices to note a few representative 
values: crio(wrci = 0-1) ~ 17, o'io(wrci = 0.3) ~ 2 and 
aio(z;rci = 0.6)«0.4 Q- 

The requirement of a non-negligible relative velocity 
then gives the decoupling time t = tdec as 

'"^^ 4.6 + (1/3) In { "^""-«^°.McV>;A 



Tdyn 



S5 



In the above, t_6 = T-dyn/lO"^ sec, and Tcoh was eval- 
uated at a terminal velocity of 0.5. The calculated de- 
coupling time is logarithmically sensitive to this choice. 



In reality the neutrons do not sharply decouple but con- 
tinue to interact with the plasma over roughly a dynami- 
cal timescale. In this sense, the t^cc appearing in Eq. (Q) 
is an "effective" decoupling time. An accurate determi- 
nation of tdoc requires solving in detail the neutron and 
proton transport equations. However, because of the ex- 
ponential decrease of density with time in the plasma 
frame, the number 4.6 appearing in Eq. (^ is only un- 
certain to approximately ±(1/3). 

Once the neutrons decouple they will have an energy 
7doc(l — Ye)M . The ratio of kinetic energy in neutrons 
to the total energy in the fireball is then 
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(Here "total" energy includes both the thermal e^/7 en- 
ergy and the bulk kinetic energy of baryons.) From this 
we see that for Yg less than 



Fe^crit = 0.46 



I ^5 
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at the time of decoupling the baryon loading problem is 
diminished. (Protons and neutrons each move with 7dcc 
at the decoupling point; thereafter, protons will possess 
larger Lorentz factors than do average neutrons). If the 
plasma remains optically thick to radiation until the en- 
ergy in radiation is converted to kinetic energy of the 
remaining protons, energy conservation gives the final 
Lorentz factor of the protons 



7 « 77 



(7) 



A simple ansatz for the condition that the fire- 
ball remains optically thick after decoupling is 77 < 
{Ye)lQ^{E/lQ^^ erg)i/3(io7 cui/R^f/^. This is obtained 
by applying the result from |^ and making the replace- 
ments ry ^ (7?/ye)(l-/n), ^ (l-/„)i; and s ^ s/^e. 
Note that even for modest values of /„, 7 can be increased 
significantly if Y^ is low. 

The above results are summarized in Fig. |l|, where we 
have plotted the smallest entropy S5 (« 1.25 r?ioo/2o) for 
which decoupling occurs as a function of Ye- For example, 
if 55 = 0.6, Tq = 2 MeV (corresponding to 77 « 100), and 
Yqo — 0.02, then the final Lorentz factor of the plasma 
after neutron decoupling (Eqs. (||), (0)) would be 7 ~ 
1500, which is 15 times larger than the standard case of 
7 = 77. As another example, consider the Ref. ||l^ values 
of T_6 = 27 and Ye = 0.1 and suppose that Tq = 10 MeV 
and S5 = 2.5 (corresponding to 77 = 2000). In this case we 
find 7 = 1.1 X lO'', an increase by a factor of 5.7. Clearly, 
the importance of this effect depends on how low Ye can 
be. 

Two conditions must be met in order to achieve a 
low Yg at the time of decoupling: (i) Yg must be low 
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FIG. 1. The smallest entropy per baryon, S5 = s/10^, for 
which decoupling occurs as a function of Ye- We have taken 

T-6O"10 = 1. 



initially and (ii) Ye must not be unacceptably raised 
during the evolution of the fireball. We can divide up 
the discussion of Y^ in this way because the initial elec- 
tron fraction depends in detail on the GRB central en- 
gine, whereas the later evolution of the fireball is gener- 
ically given by the relations in Eq. (|^). Many proposed 
GRB central engines involve neutrino heating or are 
sited in environments subject to intense neutrino fluxes 
T2I.I 
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,|18|,|19|,^,glU2|. General discussions of the 

relation between neutrino processes and the dynamics of 
outflow may be found in Refs. ||l^,^,^,^. However, 
the details of neutron decoupling are insensitive to how 
Ye is set and we are not arguing for a specific GRB site. 

The processes which have a significant effect on Ye 
in the fireball environment are lepton capture/decay in- 
volving free nucleons and inelastic nn — > np7r scattering 
(charged pion-nucleon bremsstrahlung) , 



I^e + n ^ p 
i^e + P ^ 11 + e ' 

n ^ p + + i^e 
n + n^ n + p + 7r~. 



(8a) 
(8b) 
(8c) 
(8d) 



In general, Y^ is set by the competition between the above 
processes For the range of fireball parameters of 

interest to us, free neutron decay ( |8c| ) is unimportant as 
the fraction of neutrons decaying during the evolution 
of the fireball is ~ 10~^r_6 In 77. Furthermore, as lepton 
capture is only important during the early, hot, evolution 
of the fireball and inelastic nucleon-nucleon scattering 
only occurs after neutron decoupling, the lepton capture 
and pion bremsstrahlung processes may be considered 
separately. 

In environments where neutrino heating is important 
the forward reactions ( |8a| ) and (Sb) can dominate in set- 
ting the electron fraction Integration of the rate 
equations corresponding to the lepton capture processes 



gives n/p « Xc^p/Xu.n « (Lp^(£;pJ)/(L^^(£'^J), where 
Ap^p and A^^n are the rates for the reactions in Eqs. (|8a| ) 
and ( |8b| ) , (i?p^ ) and {E^^ ) are the average energies char- 
acterizing the energy spectra of the De and neutri- 
nos, respectively, while Lp^ and L^,^ are the correspond- 
ing energy luminosities. Absent neutrino oscillations and 
flavor/type mixings, any thermal neutrino emission sce- 
nario from a compact object will yield a characteristic 
average neutrino energy heirarchy for solar mass scale ob- 
jects: {E,J « {E,J « {E,^) ^ {E,^) > (E,,) > (E,^). 
These considerations are consistent with findings in Ref. 
| ]l2| in which a hard Ve spectrum from a collapsing neu- 
tron star leads to an electron fraction in the fireball of 
Ye ^ 0.1. 

If the Ve component of the neutrino emission were to 
disappear or be greatly reduced, then the competition 
inherent in the above equations would be unbalanced in 
favor of the reaction + p ^ n + . This, in turn, 
would result in the wholesale production of neutrons. In 
fact, several schemes involving matter-enhanced active- 
sterile neutrino transformation have been proposed as a 
way of enabling r-process nucleosynthesis in neutrino- 
heated supernova ejecta: one of these involves matter- 
enhanced Ve ^ Vs and De ^ i^s 0; the other involves 
matter-enhanced conversion v^^r ^ Vs followed by an 
active-active matter-enhanced conversion f^^T- ^ i^e @. 
In either case, the intial Ve flux can be reduced by more 
than an order of magnitude and, in turn, this can trans- 
late into a substantial decrease in the initial Yg. (Just 
how low depends on central engine outflow hydrodynam- 
ics and on neutrino background effects P, |26| , p7| , p8[ .) 

If we demand that an initially low Yg not be raised 
above ye.critj consideration of lepton capture on neu- 
trons allows us to place rough constraints on the fire- 
ball and neutrino parameters. We incorporate the un- 
certainty in the initial fireball evolution by supposing 
that the relations in Eq. (||) are valid only after the fire- 
ball has a Lorentz factor 7^ and temperature Ti. Con- 
sideration of positron capture after 7 = 7i, T = Ti 

then leads to Ti < (22 MeV) (Ye.crit/'T-e)^^'^ Similar- 
ily, consideration of Pg capture on neutrons leads to 
T^^ < (40 MeV)7i (Yg^cnt/r-e)^^^ . In deriving this hmit 
we have taken the Vg spectrum to be a Fermi-Dirac black- 
body with temperature T,^^ and zero chemical potential. 
This limit could be modified or weakened if Vg flavor 
transformation occurs. 

Determining the increase in Yg due to pion produc- 
tion requires a proper treatment of neutron transport in 
the plasma. However, an upper limit on the increase is 
readily obtained by considering the extreme case where 

(i) the protons are frozen into the accelerating plasma 

(ii) non-forward n-p coUissions are assumed to result in 
maximal momentum exchange (iii) n-n coUissions are ig- 
nored except as a post-processing step to determine tt 
production and (iv) the change in Yg due to inelastic n-p 
and inelastic p-p scatterings is ignored. This simple pic- 
ture gives an upper limit on the increase in Yg because an 
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been shown to have implications for the electromag- 
netic and high-energy neutrino signatures of GRBs 
]l3| , p9|30 |. Future large volume detectors such as 
AMANDA/ICECUBE will be able to provide high en- 
ergy neutrino data on GRBs Perhaps the details 
of neutron decoupling and the associated electromag- 
netic/neutrino signature could allow a diagnostic of the 
weak interaction physics deep in GRB central engine en- 
vironments. 

It is a pleasure to acknowledge discussions with 
N. Dalai, M. Patel, X. Shi and J. R. Wilson. This re- 
search was supported in part by NSF Grant PHY98- 
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FIG. 2. The neutron distribution function /(7«) (normal- 
ized so that J fd{'yv) = 1) at several time slices as calculated 
under the assumptions given in the text. The different time 
slices correspond to a plasma 711 of 1, 1.5, 2, 2.5, 3, as mea- 
sured in a frame comoving with the plasma near the decou- 
pling point. 



exchange of any of the assumptions (i)-(iii) for more real- 
istic ones has the effect of decreasing the velocity diper- 
sion of the neutrons. By a calculation with the above 
assumptions we obtain the upper limit on the increase in 
Ye to be AYe < 10~^/yeO, where Yeo is the initial elec- 
tron fraction. Fig. (||) displays the evolution of the neu- 
tron distribution function as calculated with the above 
assumptions. The distribution function drops sharply at 
the instantaneous plasma velocity because there is no 
mechanism for boosting neutrons to higher velocity. 

The ineffectiveness of pion bremsstrahlung in increas- 
ing Ye may be attributed to the fact that the processes 
which increase Ye result from a two step process (i.e. an 
n-p scattering boosts a neutron which then inelastically 
scatters with another neutron), the decrease of a{v)v 
with increasing velocity (the product of cross section and 
relative velocity decreases by a factor of 5 as w increases 
from to c |]l^), and the late onset of this process in n-n 
scattering. Pion production does not begin until the pion 
mass threshold is reached at a relative velocity of 0.645 
and even at a relative velocity of 0.728 (center of mass 
energy 2.08 GeV), the inelastic contribution to the cross 
section is only 8% of the total Q . For AYe > Yeo our 
perturbative approach to calculating AYe breaks down 
and the increase in Ye may cause a recoupling of the pro- 
ton and neutron flows. Note that our calculation only 
gives an upper limit on the increase in Ye due to inelastic 
n-n scattering. A more careful transport calculation will 
likely show a smaller increase in Ye . It is consistent then 
to discuss decoupling at low Yg for a wide range of fireball 
parameters. The final Lorentz factor of the plasma may 
then be substantially increased for a given energy input 
and baryon load. 

Aside from the fireball energetics versus Ye issue ad- 
dressed here, neutron-proton separation recently has 
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